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Indian Standard
GUIDE FOR EQUIPMENT RELIABILITY
PART 4 PROCEDURE FOR DETERMINING AND CONFIDENCE LIMITS FROM RELIABILITY DETERMINATION

TESTING
POINT ESTIMATES EQUIPMENT TESTS

0.FOREWORD
0.1 This Indian Standard (Part 4) was adopted by the Jndian Standards Institution on 14 June 1985, after the draft finalized by the Reliability of Electronic and Electrical Components and Equipment Sectional Committee had been approved by the Electronics and Telecommunication Division Council. 0.2 This standard (Part 4) is one of the series of Indian Standards for equipment reliability testing and covers procedures for determining point estimates and confidence limits from equipment reliability determination tests. Other standards available in this series are given on page 28. 0.3 This standard is largely based on `IEC Dot: 56 (Central Office) 92 Draft, IEC Pub 605 Equipment reliability testing: Part 4 Procedures for dstermining pair t estimates and confidence limits from equipment reliability determination tests' issued by the International Electrotechnical Commission (IEC). 0.4 This standard is one of a number of Indian Standards on reliability of electronic and electrical components and equipment. A list of standards published so far in this area is given on page 28.

1. SCOPE 1.1 This standard (Part 4) gives recommended numerical and graphical methods for determining point estimates and confidence limits of reliability characteristics from equipment reliability determination tests. 3

$S : 8161 ( Part 4 ) - 4985 2. TERMINOLOGY 2.0 For the purpose of this standard, the terms and definitions in IS : 1885 (Part 39)-1979* and IS : 7690-1975.r shall apply.

as giveEa

3. LIST OF SYMBOLS
3.1 The following symbols are used in this standard. The `determination point' is the point in time or number of trials at which the point estimates or confidence limits are determined.

= location parameter in a Weibull distribution; charateristic life (or scale parameter) in a Weibul% distribution; probability density function of time to failure;

cumulative distribution of times to failure, probability of failure within time t; theoretical value of the F-distribution with v1 degrees of freedom in the numerator and 11~ degrees of freedom in the denominator at the fractile of order p; = order number failure (ti); of a failure based on the time toi

= shape parameter of a Weibull distribution; = true mean time between failures; point estimate of mean (observed value); time between failures

true mean time to failure, mean life; point estimate of mean time to value; failure (observed of

true mean time between failures for a period time (0. to);
*Electrotechnical vocabulary: items (firsf revision). Part 39 Reliability of electronic for reliability

ar,d electrical of electrxk

tMathematica1 guide to the terms arid definition; equipment and components (or park) used therein.

4
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FG (0,
13

to)

= point estimate of mean time between failures for the period of time (0, to) (observed value); = relevant number of test items or trials; = 50 percent (median) ranks of ti (see Table 1); = total number of relevant failures determination test is -based;
= true success ratio;

P,o (ti) r
R

on which

the

j,
R (t)=l-F

= point estimate of success ratio (observed value); (r) = probability of success;

t 1*m
ti

= time (or equivalent, such as distance, cycles or other quantites as appropriate); = predetermined
= relevant

test time for the mth item;

test time of magnitude order number i, recorded between the point of time the corresponding item was put on test and the point of time of failure of the item;

CP (VI

=

theoretical value of the Student t distribution with y degrees of freedom at the fractile of order p;

T*

= accumulated relevant test time up to the determination point when this does not coincide with a failure; = relevant test time accumulated for all items up to the determination test point when this coincides with a failure; = true constant failure rate; = point estimate of failure rate (observed value);
=

true standard deviation in normal distribution;

= point estimate of the standard deviation; and = theoretical value of the X" distribution with v degrees freedom at the fractile of order P.
5
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TABLE 1 MEDIAN RANKS ( P 50 IN PERCENT ( Clauses 3.1, 6.3.1, 7.1.1, 7.1.2 and 7.2.2.1 ) SAMPLE SIZE ( n ) ~~_~~-_~_________h______~_~~__~~~~~ 5d.O 2922 70 7 20 36 50.0 70.3 15.9 4 38.5 61.4 84.0 12.9 5 31.3 50.0 68.6 87.0 10.9 6 26.6 42.1 57.8 73.5 89.0 974 22.8 34.4 50.0 63.5 77.1 90 5 )

ORDER No?)

6 7

8

8.3 8 20.1 32.0 44 0 55.9 67.9 79.8 91 7

1:
ORDER No. (I)

9. I:.; 28.6 39.3 50.0 60.6 71.3 82.0 92-5

10 6.6 16.2 25.8 35.5 45.1 51.8 64.4 74.1 83.7 93.3

SAMPLESIZE ( n )

~-------____-_--11 6.1 14.7 23.5 32.3 41.1 50.0 58'8 67.6 76.4 85.2 93.8 12 5.6 13 5 21.6 29.7 37 s 45.9 54.0 62.1 70.2 78.3 86.4 94.3

h_______--_-_-_-_-_-_?

: i
5 6 7 8

1;: :a 13 14 15 16 :; :`o

20.0 27.5 35.0 42 5 50 0 57.4 64.9 72.4 79.9 87.4 94.8

::1 12-5

14 4-S 11.7 IS 6 25.6 32.5 39.5 46 5 53.4 60.4 67 4 74.3 81.3 88.2 95.1

10 9 17.4 23.9 30 4 36.9 c3.4 50.0 56.5 63.0 69.5 16.0 82 5 89 0 95.4

25

16 . 140.5 16.3 22.4 25.5 34.7 40.8 46.9 53.0 59.1 65 2 71.4 77.5 83.6 89.7 95.7

18

9.6 154 21.1 26.9 32 7 38.4 44.2 50.0 55.7 61'5 67.2 73.0 78'8 S4.4 90 3 96.0

14.5 20 0 25.4 30.9 36.3 41.8 47.2 52.7 58.1 63.6 69 0 74.5 79.9 85'4 90 s 96.0

?.I

8.6 13.5 1x.9 24.1 29.3 34.4 39 6 44.8 50.0 55.1 60.3 65.5 70.6 75.8 81.0 86.1 91.3 96'4

NOTE -

For higher values of IZand i use formula: i - 0.3 P," = 100 n

3.4 8.' 13.1 IS.0 22.9 27.8 32.7 37.7 42 6 47.5 52.4 57.3 62.2 67.3 72.1 77 0 81 9 S6.S 91.7 96.5

26

4. APPLICABILITY 4.1 The methods of calculation in this standard nlay be applied at any time or after any number of trials of the reliability testing. The more information .that is available when the estimates and limits are determined, the higher is the precision obtained. Existing data from earlier tests or field observations may be acceptable, provided the data are sufficiently complete, well establised and applicable to the situation. The standard is not applicable to the pooling of data from different sources, for example, original data from tests under different conditions. 6

IS : 8161 ( Part 4.2 Wherever `time' replaced by distance,

4 ) - 1985 be

is used in this standard, this variable may cycles or other quanties as may be appropriate.

4.3 If a constant failure ra!e, a mean time to/between failures, or a mean failure rate for a given period of time is to be estimated from the reliability determination test, the methods in 6 or 7 should be used. These methods are applicable for tests based on time. The method to be used depends on what assumption can be made regarding the failure rate depe.ldence on time and the statistical distribution of time to/ betwee; failures. 4.4 In case of an assumption of a constant failure rate characterising an exponential distribution of times to/between failures, the method given The test items may be repaired and put back on in 6 is to be used. test or not repaired. The graphical method (6.3) is, however, not applicable after the First failure of repaired test items. In case of an assumption of a non-coi;stant failure rate described by a Weibull distribution or a normal (Gaussian) distribution of times to failure, the methods given in 7 should be used. These methods are not applicable after the first failure of repaired test items. 4.5 If the success ratio is to be estimated from a reliability determination test, the methods in 8 are to be used. These methods are based on the number of test items or trials, which are classified as either failed or no:l-failed. The method to be used depends on the number of failures. If the number of failures is greater than 9 the normal distribution may be used, in other cases the binominal distribution shall be used. 5. FEATURES LIMITS OF THE POINT ESTIlMATES AND CONFIDENCE value to represent the such as the failure rate. is the `observed' value : 7690-1975`: 1.

5.1 A point estimate is a single numerical unknown true value of a statistical parameter, Generally, the point estimate considered here [ as defined in IS : 1885 ( Part 39 )-1979" and IS

5.2 The confidence limits define a confidence interval around the point estimate where the interval includes the true value of the parameter being estimated with a certain probability, the confidence level. 5.3 The confidence interval available. The information and the number of failures will be narrower when more information is used is the accumulated relevant test time for time based tests, and the number of test

*Electrotechnical vocabulary: Part 39 Reliability of electrcn'c and electrical items (first uevisio~z ). tklathematicnl guide to the terms and definitions for reliability of electronic equipment and components (or parts) used therein. 7
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or trials

and

the number of failures for tests based on items or

.trials. 5.4 The confidence interval may be one-sided or two-sided. In case of .a one-sided confidence interval, a lower or an upper confidence limit is .given for the characteristic. In case of a two-sided confidence interval, both a lo\\,er and an upper limit are given. 5.5 The preferred confidence level is here chosen to be 90 percent. The confidence intervals according to this recommendation therfore include the true values of the characteristics with 90 percent probability. If other confidence levels are desired, the formulae in this standard may be used with appropriate values of the test statistics taken from statistical tables. 6. CONSTANT FAILURE RATE

6.0 If the true failure rate is assumed to be constant, the numerical and graphical methods in 6.1, 6.2 and 6.3 may be used for estimating failure rate and mean time between failures (for repaired items) or mean time to failure (for non-repaired items). The methods are applicable both .to repair and non-repair situations. These estimates are based on the total number of relevant -failures and the accumulated relevant test time UP to the determination point: 6.0.1 The validity of a constant failure rate assumption should prefer.ably be tested in accordarce with IS : 8161 (Part 6)-1983* before point .estimates and confidence limits are calculated. For application For application to mean time between failures, M.T.B.F. to mean time to ( first > failure,
nlF = mF

= 111.

or M.T.T.F.:

M.T.T.F.

= m. -L

For application

to failure rate h =

$6.1 Time Terminated Tests 6.1.0 The accumulated relevant test time up to the determination with Appendix A. The point, T*, may be determir,ed in accordance ,same formulae are used for test with and without replacement or repair. The total number of relevant failure, r, up to the determination point ;is counted.
*Guide for equipment reliability testing : Part 6 Tests for validity of a constant

:failrlre rate assumption.

8
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5.1.1 Point Estimates

The point estimate (observed value) of the failure rate is: i = T< The point estimate (observed value) of the mean time to or between failures is: T" h
in = -;

If no failures are observed up to the termination point, that is, r = 0, the following point estimate of the failure rate is recommended: ^h = 1/3T* Too much reliance should not be placed zero or any small number of failures. on estimates based upon

Care should be taken in planning the tests such that T* will be sufficiently large in relation to the expected mean time to failure. NOTE - This recommendation is based upon an investigation by E.L. Welker
and h$. Lipow `Estimating the exponential events' in Proceedings of 1974 Annual posium, pages 420-427. failure rate from data with no failure Reliability and Maintainability Sym-

6.1.2 ConjLlence Interval.9
(j.1.2.1 The confidence limits for the true failure rate, h, with the confidence level of 90 percent are shown below (see Table 2 for values of the X" distribution).

a) One-sided confidence
X2

interval, upper limit:
2)or h ( Xz~.9~2( Z&f 2)

X</?b) Two-sided i or

0.~0(2r+ 2f

confidence
2r

interval: <AA uh2.95 ( 2r + 2)
2r

x20*05 (2v)

X'cl.,; ( 2r) 2T*

< A < X" 0.95 ( 2r + 2) 22"" limit with

If no failures are observed only the one-sided confidence an upper limit can be defined.

IS : 8161 ( Part 4 ) - 1985
TABLE 2 CM-SQUARE DEGREESOF FREEDOMY 2 4 6 8 :; 14 16 :"o 22 24 26 28 30 32 34 36 :: Xeo.lJ5(v 1
.103

DISTRIBUTION XZrJ.,,( y ) 4,605 7.779 10.65 13.36 15.98 18.55 21.06 23.54 25.99 28.41 30.81 33.2 35.56 37.92 40.26 42.57 4488 47.19 49.5 51.81 54.08 63.17 65.42 14.4 76-63 85.53 a7 74 96.58 98.78 107.57 109.76 118.5 120.65 129.25 131.4 140 142.15 226.02 +1.29

FRACTILES ) Xa,.,J V ) 5.991 9.488 12.59 15.51 la.31 21.03 23.69 26.3 28.87 31.41 33.92 36.42 38.89 41.34 43.77 46.17 48.57 SO.96 53.36 55.76 58.11 67.51 69.82 79.08 81.37 90.53 92.8 101.88 104.13 111.15 115.39 124.34 126'53 135.3 137.5 146.27 148.46 233.99 t1.64 accurate,

( Clauses 6.1.2.1, 6.1.2.2, 6.2.2.1 and6.2.2.2

' 42
50

52

6i.i
62 70 72 80
90

si

92 100 102 110 112 120 122 200
Z=

*711 1.635 2.733 3.94 5.226 6.571 7.962 9.19 _IO.85 12.34 13'85 15.38 16'92 18.49 2009 21.7 23'3 24.91 26.51 28.16 34.76 36.45 43.19 44.9 5174 53.47 60.39 62.14 69 13 70 89 77.93 7974 86.96 88.77 96 97.81 168.28 -1.64

NOTE 1 - Linear interpolation of intermediate values is sufficiently Values for v=2r+2 for various integer r values are included.

NOTE 2 -

For higher

values, use: percentage of the standard

Xp ( Y )=[ (z+2/2v-1 )"j!2 where z is the corresponding normal distribution given at the foot of each column.

10

IS : S-161 ( Part 4 ) - 1985 6.1.2.2 The confidence limits for the true mmeantime to or between failures m, with a confidence level of 90 percent are shown below (see Table 2 for values of the X" distribution).

a) One-sided confidence interval, lower limit:
m> GI

2r
X`o.90 (

2r + 2)

Or m '

%`)O.QO 9;: + 2)

b) Two-sided confidence interval:

G
or

x2

tmtcl 0.95 6-r + 2 1,

X2o.05 Fr2r 1

XL 0.95 $"+

2)

<m<

2T*
x"0.0~ ( 2r >

If no failures are observed, only the lower limit can be defined. The 90 percent confidence limits as a, function of the number of failures, r, are shown in Fig. 1 and Table 3. The limits are here expressed by the appropriate multiplier times the relevent point estimate x^and rn?
NOTE- The multipliers may be used in the planning of life tests by finding out the approximate number of failures for a given precision of the estimate. The number of failures is required accumulated test time, T*? for the required approximately T*=y/A, where the failure rate A has to be assumed from earlier experience.

6.2 Failure Terminated Tests - The accumulated relevent test time up to the determination point, T,, is determined in accordance with Appendix A. The same formulae are used for tests with and without replacement or repair. 6.2.1 Point Estimntes failure rate is:
-

The point ~=_.$-

estimate ( observed

value > of the

The point estimate ( observed value ) of the mean time to or between failures is:
Tr

Too much reliance should not be placed on estimates based upon a small number of failures. Care should be taken in planning the -tests such that T* will be sufficiently large in relation to the expected mean time to failure. 11
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r4ULTIPLIER 1 FOR FAILURE RATE

ONE-SIDED

UPPER

LIMIT

'

WO,-SIDED

INTERVAL

MULTtPLlER 1

FOR

FAILURE

RATE

TWO-SIDED

INTERVAL

I ONE-SIDED LOWER LIMIT

FIG.

1

THE 90 PERCENTCONFIDENCE LIMITS AS A FUNCTION OF NUMBER OF FAILURES FOR TIME-TERMINATED TESTS

,12

IS : 8161 (apart 4 ) - 1985
TABLE 3 MULTIPLIERS TO OBTAIN CONFIDENCE LIMITS

( Clause 6.1.2.2 ) MULTIPLIERS FOR MTTF/MTBF MULTIPLIERS FORFAILURE RATE h_______~ h~~~~~~~~ c------r------R

percent )
:: 3 t 6 7 9" 10 :1 Cl3 14 E 17 1s 19 20 :; 35 40 45 50 55 60 19.417 5.625 3.669 2.927 2.518 2.296 2.13 2.009 I.916 1.843 1.782 1.69 1.732 I ,654 1.622 1.592 I 566 1.545 1,525 1.508 1.438 1.389 1.352 1.324 I .301 1.283 1.264 1,249

M(5

M ( 90

M ( 95

M(5

M ( 90

percent )
-257 -375 *449 *5 *539 .569 .594 .615 .633 *649 .662 -685 .674 .695 .704 $3 .727 .733 `739 .764 `782 .797 *so9 .819 *828 .837 `844

percent )
.21 *317 -386 .436 *475 `506 .532 a554 `573 .589 *604 .628 .617 .639 .649 .658 `667 -674 *681 .688 *716 .I37 -754 .768 `779 .79 ZOS

percent )
*051 a177 -272 *?41 .394 .435 -469 .497 -521 .542 .56 .591 -577 .604 `616 .627 .638 .64.7 .655 .662 ah95 .719 -739 *754 `768 *779 :`8"

percent )
3.889 2 662 2.226 1.997 1.854 1.755 1,681 1,624 I .578 1.54 1'509 1.458 I .481 I.437 I.419 I.402 1.387 I.375 I.363 1,351 I.308 I 257 1.253 1.234 1.219 1.206 I.194 1.184

M ( 95 percmd ) 4743 3.147 2.584 2,288 2.103' 1.979 1'878 1.804 1.745, 1 ,G95 1.655 :::f I.563 1.53% I.516 3.498. 1.482 1.467 1.452 I .396' 1.356. 1.325 1.301 I 282 : F I .337

5.2.2

Confidence

Intervuls

6.2.2.1 The Confidence confidence level of 90 percent of the X" distribution ):

limits for are shown interval,

the true failure rate with a below ( See Table 2 for vnlu@s. upper limit:

a)

One-sided

confidence

b)

Two-sided

confidence
2r

intxval: XJ."O.$IS ( 2u )
2Tr

;ALo'i or

( 2l. ) < x <;

x20.,5 ( 2r 1 < x < 2T1 13

%.s~ ( 2r 1
21r

SS : 8161 ( Part 4 ) - 1985 6.2.2.2 The confidence limits for the true mean time to or between failures with a con5dence level of 90 percent are shown below, ( see `Table 2 for values of the X' distribution ). a) One-sided confidence interval, Orm ' interval:
< 111 <

lower limit: 2Tr X'O.~O ( 2r )

,S+&) Two-sided
h

2r
0.20 ( 2r ) confidence 2r X.`)o.1)5 ( 2r )

IT1

2r nT Z'o.oj ( 2r

j

or 2Tr
X?O.~~ (

2Tr
2r )< m < Xzo.05 ( 2r j

$.3

Graphical

Method method is presented for *estimating failure rate and mean failure or mean time between failures, based on the semiThe tests from which the observations are taken need not all test items have failed. be used only for times to first failure The method will give point estimates of deviations from a constant failure observed. For

63.0 A `time to first Iog paper. go on until

The semi-log paper shou!d ,and for not less than 4 failures. .and may also give an indication rate. 6.3.1 A number of items, brach failed item, thz relevant The times ti are ordered

II, are tested and r failures test time tl is recorded. in magnitude:
tr

tl < i`z < t:... . ..<

cd9 a semi-log pnpcr ( see Fig. 2 ), the ti values are plotted along .+& linear scale and the reciprocal of one minus the fractional median rrrnk, P,5o ( ti)/lo0 along the logarithmic scale. Values of median ranks may be found from Table 1. If a constant failure rate assumption holds, the plotted points will fit v,ell to a straight line passing through the the line, the central points shou!d point, fi=O, ratio = 1. In drawing dominate the determiniration of the slope. This being so. then the estimate of the mean time to fai!ure equals the t valise on the time axis that corresponds with the ratio 2'72 on the vertical axis. The estimate `of the failure rate is then the inverse of the t va!ue.

14
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REI."V&NT

TEST .TIME FOR EACH ITEM

FIG.

?

EXAMPLE

OF USE OFSEMI-LOG PAPER
)/lUO

1 1 - P;o
(ti

;;l=-;;- 1 A 3f the plotted points cannot be Ctted by a straight line, the failure I!1 such a case, other n:ethods as described in sate may not be constant. 7 may apply. 6.3.2 Example 6.3.2.B Figure 2 has been drawn with data col!ected by putting 10 items in reliahiiity test which was terminated when 8th failure was observed. The observation and calculations are given below: l-&O C ti) hours ) 100 29 6'6 1 1 07 29.9 16'2 2 1'19 30'6 25 8 3 1'38 4 32'4 35.5 1'55 33'0 45 1 1'82 :, 35'3 54'8 2'21 7 36'1 64'4 2 81 74'1 40' 1 8 3'86 6.3.2.2 The above dnta is plotted on semi-log paper in Fig. 2 as described in 6.3,P and observed value for estimated failure rate at the
(

Failure O#"u%r

Test Time (ti) and ith Fnilure

Pso (li)"

1

percent

ratio value of 2'72 is equal to &
*Obtained from Table 5.

h = 0`027 4 failures/hour. 15

IS : 8161 ( Part 4) - 1985 7. NON-CONSTANT FAILURE RATE

7.0 If the assumption of a constant failure rate is not valid and the distribution of times to failure follow a Weibull or a normal distribution, the methods given below are applicable. The validity of any of these distribution assumptions should preferably be tested in accordance with IS : 8161 ( Part 8 )* or using the probability papers described in 7.1 and 7.2. 7.Q.l The life history of the test items shall be similar when they are put on test. The methods of evaluation are applicable only to the times to first failure of each item. The test does not, however, need to be continued until all items have failed, except for the method in 7.2.1.1. 7.0.2 For each failed test item, the relevant test time, tr, is determined. The time ti is defined as the relevant test time between the point of time when the item i was put on test, t = 0, and the point of failure, t = tiThe times ti are ordered in magnitude
tl < t:! -=z ta < . . . . . . < t,

Using these times and a Weibull probability paper or a normal probability paper it is possible to estimate the parameters of the distribution and at the same time obtain an indication of how well the times are fitted to the assumed distribution. 7.1 W&bull Distribution 7.1.0 The cumulative distribution is: F(t)=l-Jt;%; > a, b > o, k > o function for the Weibull distribution

7.1.1 A graphical method for point estimates is presented below. In this method, the Weibuil probability paper is used. This paper is so designed that the cumulative distribution function of times to failure which follow a Weibull distribution, will become a straight line. See Fig. 3 which is an example of such a paper and how it can be used. If the time scale on the abscissa is not directly applicable, transformed by a scale factor lo", where c is an interger. it can be

The observed ti values are plotted at the corresponding median rank values ( P,o) ( tr ) on the F ( t ) scale. Median rank values can be found in Table 1.
*Guide for equipment reliability testing: Part constant failure rate assumption ( urzrler consideration 8 Tests
).

for validity

of

anon-
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RELEVANT FIG.

TEST TIME FOR EACH ITEM (t 14

3

EXAMPLE-OFUSE OF WEIBULLPROBABILITY:PAPER

7.1.2 test is plotted for n =

The example of Fig. 3 assumes a terminated after r = 8 failures. points and in the table of median 10 is used up to and including i =

test sample n = 10 when the Therefore, there are only 8 ranks ( Table 1 > the column 8 only.

7.1.3 If a Weibull distribution with a = 0 holds then the plotted points can be fitted will by a straight line. If the points follow a concave curve seen from the abscissa, the location parameter a has a positive value. If the curve is conves then a is negative. For a curve the following procedure may be used:

Step 1: Plot curved data as described above and illustrated in Fig. 3. Step 2: Draw two parallel horizontal lines through the extreme failure points, and a third horizontal parallel line at the linear mid-point. Step 3: Draw three lines perpendicular `to the t-axis from the points of intersection of the three horizontal lines and the curve. St.ep 4: Identify the lowest, and ta respectively. middle, and highest t value as tl, tm

17
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Step 5: Use the following ( th L+_
( th -

formula to calculate:
tm ) ( tm ~tm ) ( tm tl tl > >

Step 6: Subtract ^a(including sigh of a^) from each tr value. Step 7: Replot the data; the points should now approximate a straight line, a shifted Weibull distribution of failures. 7.1.3.1 The data plotted in Fig. 3 was collected on 10 items kept on reliability test and the test was terminated at 8th failure. The observed failures and calculations as described above are given below: Failure Order : : 5 6 s' Test Time, at ith Failure ( hours) 400 510 620 720 820 920 1 050 1 300
th =

Medium Rank Value percent* 6'6 16`2 25'8 35'5 45'1 54.8 64.4 74' 1 1 300 and
tm =

7.1.3.2 From above data, tl = 400, hence a = 294.8 ( m 300 ).

620 and

7.1.4 TO estimate the shape parameter, k, a line is drawn through the centre of the circular scale parallel to the fitted straight line. The value of k is then read off the circular scale at the intersection. The estimate of the characteristic life or scale parameter, b, is equal to the t value of the intersection of the fitted straight lice and dashed line on the vertical axis at 6396. 7.1.5 To estimate mean life mF refer k values the value of F ( t ) when t ordinate scale of F ( t ) in Fig. 3, mF plot ( not the straight line derived The original plot may also be used the population having times to failure to Fig. 4. This gives for different Then by reference to the can be read off from the original from a ) at the abscissa scale. to estimate the fraction F ( t ), of less than a given t value.
= mF.

7.1.-6 Other available types of graph paper, while basically similar, may use different methods to estimate the parameters asinstructed.
*Taken from Table 5.

18
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7.2 Normal Distribution 7.2.0

The probability

density function of the normal distribution

is:

This distribution 7.2.1 Numerical

is useful for times to failure,
Method

provided

that

mF

>

30.

7.2.1.1. Point estimates - The numerical estimates of the mean time to failure and standard deviation presented here are applicable only when all the 12test items are tested until they fail. 19
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a) The point estimate of the mean time to failure is:
h

mF =n

1

c
i=l

n

Ii

b) The point estimate of the standard deviation is:

I
A=
0

d

i=l

tt

( ti -

m^F >" =

I

n-l

d

n A i=l

fi3 - &ti 1" `I7 n ( Iz 1)

7.2.1.2 Confidence intervals - The confidence limits for the true mean time to failure, iTiF, with a confider.ce level of 90 percent are shobtn below [ see Table 4 for values of the Student's distribution tl, ( 1:)l.

TABLE

4 STUDENT

t DISTRIBUTION
DEGREESOF FREEDOM v I tLv!lo (VI

D~ciRtEs OF

tcl.90

to.95

FREEDOM Y
2 3 4 5 6 7 8 9 10 12

Cv)
1.89 1.64 1.53 1.45 1'44 1 41 1.40 1.38 1.37 1.36

(Vj
2.92 2.35 2.13 2.02 1.94 1.89 1.86 1.83 1.81 1.78 ~

14 16 18 20 25 30 4o 60 100 cc

1.34 1.34 1.33 1.33 1.32 1.31 1.30 1.30 1.29 1'28

1.76 1 75 1.73 1.72 1.71 1.70 1.68 1.67 1.66 1,64

KOE -

Linear

interpolation

for intermediate

Yvalues

is sufficiently

accurate.

20
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a) One-sided confidence interval, ?wer h
??lF > mF tO.$O(n 1 ) fn

limit:

b) Two-sided
niF _

confidence interval: h
to.gs ( n 1 > ,yn. <
mF < $?F +t0.B5(n1)"

h
An

`7.2.2 Graphical Method

RELEVANT

TEST

TlME_FOR

EA-CH

TIME

In example, n= 10, r= 8 Failures occurred times in hours as plotted mF= 1 500 FIG. 5 EXAMPLE OF USE OFNORMALPROBABILITYPAPER
21

IS : 8161 ( Part 4 ) - 1985 7.2.2.0 The graphical method is based on the normal probability paper. This paper is so designed that the cumulative distribution function of times to failure, which follow a normal distribution, will become a. straight line. see Fig. 5. It is applicable also in the case that not all items have failed. 7.2.2.1 The observed ?I values are plotted at the corresponding median rank values, pjo ( ti ), on the F( t ) scale. Median rank va!ues t - Q'3 where ___ ~ 12 -I- 0'4 n is the number of test items. The example of Fig. 5 assumes a test sample n=lO with the test terminated after Y = 8 failures. In such a case, where the test is-terminated becore all items have failed, only the median ranks up to and including i = P are plotted. In drawing the line the points between the probabalities 20 percent and 80 percent should principally determine the line. may be found in Table 1 or approximated by the ratio 7.2.2.2 The estimate of the mean time to failure, the t value of the point on the straight line, which equal to SOpercelnt.
~ZF, is equal to has an ordiliate

7.2.2.3 The estimate of the standard deviation, a^, is equal to the difference between the t value of the point on the straight line which has an ordinate equal to 84 percent and the estimated mean time to, failure. 7.2.2.4 The fitted straight line may also be used to estimate fraction of the population below or above a given time to failure. 8. SUCCESS 8.0 The required RATIO ratio is the probability that or that a trial will be successful an item will perform a under stated conditions. the-

success function

When testing for success ratio, the test items or trials are to be classified as either failed or non-failed. Reusable devices may be repaired between successive trails, provided that the state and performance are the same at the start of all trials. The point estimate ( observed success ratio ) and the confidence interval for the success ration are independent of how the test is terminated at a fixed number of trials, a fixed number of successes or a fixed number of failures. 8.1 Point Estimate - The point estimate of success ratio is equal to the num.ber of successes divided by the total number of trials or test items:
esn-r n

8.2 Confidence

Intervals for the true success based on the binomial 22 ratio R with a distribution,. are'

8.2.0 The confidence limits confiderAce level of 90 percent,

IS : 8161 ( Part 4 ) - 1985 shown below. ( See Tables 0'9 and 0'95 fractiles. )
TABLE V,\Vl 2 4 6 8 i; f ;t 18 20 30 2: 120 co
NOTE

5 and 6 for Values of the Fdistribution
20 0.9 FRACTlLES 30 40 9.47 3.80 2.78 2.36 2.13 1.99 1.89 1.81 1.75 1 71 1.57 1.51 1.44 :`:o' 60 9.47 3.79 2.76 2.34 2.11 1.96 1.86 1 78 1.72 1.68 1.54 1 47 I .4O 1:;; 120 9.48 3.78 2.74 2.32 2.08 1.93 1.83 1.75 1.69 1 64 1.50 1.42 1.35 1.26 1.17

for

5 F DISTRIBUTION 6 9.33 4.01 3.05 2 67 2.46 2.33 2 23 i.1: 8 10

* 2 9.00 4.32 3.46 3.11 2.92 2 81 2.73 2.67 2.62 2.59 2.49 2.44 2.39 2.35 2.30 - Linear

4 9.24 4.1 I 3.18 2 81 2.61 2.48 2.39 2 33 2.29 2 25 zb:

co 9.49 3 96 2-72 2 29 2.06 1.90 1.80 1.72 1.66 1.61 1.46 1.38 1.29 1.19 1.00

2 09 1'9s 1.93 2.04 1.87 1.99 1.82 1'94 1.77 interpolation

9.39 9.44 9.46 9.37 3.92 3.84 3.82 3.95 2.91 2.84 2.8!) 2.98 2.54 2.42 2.38 2.59 2.32 2.20 2.16 2.38 2.19 2.06 2.01 2.24 2 10 1.96 1.91 ?.I_5 2.03 1.89 1.84 :Yi" I.98 1.84 1.78 0; l-94 1.79 1.74 2 1.82 1.67 1.61 1.88 1.76 1.61 1.54 1.83 1.71 1.54 I .48 1.77 1.65 1.48 ;4fi 1'72 1.60 1,42 1.67 for intermediate v values is

.

j1 L

sufficiently

accurate.

TABLE v,\v, 2 4 fi 8 10 12 14 16 :: 30 40 1% 00 2 19.00 6.94 5.14 4.46 4.10 3'89 3,74 3.63 3.55 3.49 3.32 3.23 3.15 3.07 3.00 4 19.20 6.39 4.53 3.84 3.48 3.26 3.11 3,Ol 2.87 2.93 2.69 2.61 2.53 2.45 2.37 6

6 F DISTRIBUTION 8 19 40 6.04 4 15 3.44 3.07 2-85 2.70 2.59 2 2.51 45 2 27 2.18 2.10 2.02 1.94 10 19.40 5.96 4 06 3 35 2.98 2.75 2.60 `2'%? 2.35 2.16 208 1.91 I.99 1.83 20

-

0'95 FRACTILES 30 19.50 5.75 3.81 3.03 2.70 2.47 2.31 2'19 2.11 :::r: 1.74 1.41 1.65 1.46 40 19.50 5 72 3.77 3.04 2.66 2.43 2.27 2.15 2.06 1.g9 1.79 1.69 1.37 1.59 1.39 60 13.50 5.69 3.74 3.01 2.62 2.38 2.22 2.11 2.02 1.95 1.74 1.64 1.32 1.53 1'32 120 19.50 5.66 3.70 2.97 2.58 2.34 2 18 2.06 1.9') 1.97 1 68 1.58 1.47 1.26 1.22 co 19.50 5.63 3.67 2.93 2 54 2-30 2.13 2.01 I.92 1.84 1.62 1.51 1.39 1.19 1'00

19.30 6.16 4.28 3.58 3-22 3.00 2 85 2.74 2.66 2.60 2.42 2.34 2 25 2.18 2.10

19.40 5.80 3 87 3.15 2.77 2.54 2.39 2'28 2.19 2.12 1.93 1.84 1.66 1.75 1.57

NOTE -

Linear interpolation

for intermediate

v values is sufficiently

accurate.

8.2.1

One-sided
R>

confidence
12 r +

interval,
n -r
( r +

lower

limit:

1 ) FgO ( rl, ~2 >

with

vl=2(r+l)andvz=2(n--) 23
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I? -

confidence interals:
r

n -

r +

( r -i 1 ) FS5 ( Y,, yg ) < ' <

( n - r -!- I 1,FQ, ( vl, vs ) 1 + ( II - r -t 1 ) FD5( vlr yyt)

with

v1 = 2 ( r + 1 ) and Ye = 2 ( IZ - r >at lower limits, and - r -I- 1 ) and Y.~ v,=2(n = 2r at upper limits

8.3 Use of Nomographs 8.3.1 Nomographs are presented for determining one-sided and twosided confidence intervals of success ratio, see Fig. 6 and 7. The nomogranhs are based on the total number of test items or trials, n, and the total ;lurnber of observed failures, r.

T

r = NUMBER OF FAILURES

n = TEAL

NUL(BE~J c OF TEST ITEMS OR TRIALS

r
0~0001

'

0.001

o-01

04

FAILURE

RATIO,

r/n

_C

FIG. 6 8.3.2 Enter the applicable nomograph on the Go vertically up to the curves labelled with failures, r. Then go horizontally to the ordinate failure ratio limits, Q. The confidence limits of complements, that is, R = 1 - Q.

abscissa at the r/lz value. the observed number of and to read the relevant the success ratio are the

24
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0.001

0.01

FAILURE RATIO,r/n

V

FIG.7
`1lIsr
ITEM `+LJWSER

r

FAILURE

RELEVANI TIMES

TEST

1

I

FAlLURi NUMBER

I
h CALENDER
TIME

k-1

05

FIG.

8
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APPENDIX
DETERMINATION

A

( Clauses 6.1.0 and 6.2 )
OF ACCUMULATED RELEVANT TEST TINIE

A-l. The relevant test time of each individual test item can be measured by elapsed time indicators on the test item. In this case the accumulated relevant test time 3k at the kth failure is the sum of indicator readings.
n Tr~=c where 12 = the total
tk m =

rn =

tk,m 1

number

of test items,

and m up a decision to the point kth not

the indicated relevant test item of number failure in the total number of test items. accumulated relevant with a failure is: T* = ;
m=l

A-2. The coinciding

test

time

T*

at

tm*

where G = the indicated point. test time of item number
m up to the

decision

A-3. If the relevant test time is recorded by other means, the accumulated relevant test time TI,, at the kth failure should be calculated by the This includes the accumulated relevant following iterative formula. test time up to the failure k - 1 and the relevant test times elapsed in the interval between the failures k - 1 and k. Tk = where n = the total number of test items, and tm,j = the jth period of relevant test time of test item number m after the failure k -- 1 in the total number of test items. A-3. The intermissions may be caused technical or administrative reasons. may vary from item to item. by failure k - 1 or by any other The number, j, of intermissions

Tic - 1 +

l ZJ, tm, j WI= 1

26

IS ~: 8161 ( Part 4 ) - 1985 A-4. Figure A-5. The coinciding 8 explains the numbering of periods of relevant test time. point not

accumulated relevant test time with a failure is in this case:
T* =

T* at a decision

Tr +
m

i
==

zj
1

tm j
,

where relevant test time up to the latest T, = the accumulated before the decision point, and t III, j failure

= the jth period of relevant test time item number 771after the latest failure in the total number of test items.

A-6. The above formulae are applicable also in case the test items are con-repaired equipmects, except that relevant test time does not exist after the first failure of each test item.
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